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Abstract. This paper presents an analysis of the correla- 
tion between the far-infrared (FIR) and the radio emission 
of starburst galaxies. Data for interacting galaxies, many 
of which are undergoing a starburst, and for normal galax- 
ies have been analysed and compared in order to test for 
any influence of the star-formation activity on the ratio be- 
tween the FIR and the radio emission at 2.4 GHz, <?2.4GHz- 
No statistically significant indication for such an influence 
has been found: There is neither a significant difference 
between g2.4GHz of the two samples nor a dependence of 
this ratio on the starburst strength. This observational 
fact is unexpected because of the different physical condi- 
tions and the short time-scale of the star-formation activ- 
ity in a starburst. In order to interprete the observations 
the FIR and the radio emission during a starburst have 
been modeled. The following conclusions could be drawn 
from the observed constancy of <?2.4GHz: A strong and fast 
(« 10 7 yr) increase of the magnetic field at the beginning 
of the starburst is required in order to maintain a constant 
<72.4GHz- Otherwise the strong Inverse Compton losses that 
are due to the intense radiation field in a starburst would 
lower the synchrotron emission drastically resulting in a 
value of (?2.4GHz significantly higher than the observed one. 
Furthermore, the time-scale of the variation of the star- 
formation rate has to be longer than some 10 7 yr. Fol- 
lower values the different time-scales of the FIR and the 
radio emission produce large fluctuations of (72.4GHz ■ 

Key words: Galaxies: starburst - galaxies: magnetic 
fields - infrared: galaxies - radio continuum: galaxies 

1. Introduction 

Starburst galaxies are found to follow the tight correla- 
tion between the far-infrared (FIR) and the radio emis- 
sion (Wunderlich et al. 1987; van den Driel et al. 1991; 



Franceschini et al. 1988a; Condon et al. 1991). Here, the 
correlation is even more unexpected than in normal galax- 
ies (see Lisenfeld et al. 1995 for a discussion of normal 
galaxies) because of the following reasons. In starburst 
galaxies the different time-scales involved in the produc- 
tion of the radio and the FIR emissions become impor- 
tant: Whereas the FIR emission and the thermal radio 
emission are instantaneously related to the stellar radia- 
tion, the synchrotron emission shows a time delay because 
it takes some 10 6 yr until the first supernovae (SNe) oc- 
cur, and some 10 7 yr till the maximum of SN activity is 
reached. These time-scales are comparable to the duration 
of starbursts of about 10 7 — 10 s yr (Rieke et al. 1980; van 
den Broek et al. 1991; Bernlohr 1993) and should there- 
fore result in a time delay between the maximum emis- 
sion of the FIR and the radio emission. Furthermore, dur- 
ing a starburst the intensity of the radiation field changes 
with time in a significant way. Little is known about vari- 
ations of the magnetic field energy density. The ratio of 
these two parameters that determine the major radiative 
losses of the cosmic ray (CR) electrons (inverse Compton 
and synchrotron losses) is crucial for the total synchrotron 
emission of a galaxy. 
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The aim of this paper is to understand why, in spite of 
the above considerations, the FIR/radio ratio is observed 
to be very constant in starburst galaxies, and, more im- 
portantly, to use this observational fact to find constraints 
for starburst models. The paper consists of two parts: In 
the first part (Sects. 2 and 3), we present data for two 
large samples of galaxies, a starburst sample and a con- 
trol sample, in order to investigate whether the FIR/radio 
ratio shows any dependence on the star-formation activ- 
ity. In the second part (Sects. 4-6) we present a model for 
the FIR and the radio emission in a starburst in order to 
interprete the observational results. Section 7 summarizes 
our conclusions. 
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2. The sample 

Both the sample of starburst galaxies and the control sam- 
ple are selected from the list of UGC galaxies which are 
surveyed with the Arecibo 300m telescope at 2.4 GHz 
(Dressel & Condon 1978) and were covered by IRAS in 
FIR . These galaxies form a complete sample (down to 
me = 14.5) of a local galaxy population, which have been 
studied by Franceschini et al. (1988a, b) for the radio and 
FIR luminosity functions. 

It has been well established that non-Seyfert Markar- 
ian galaxies (Huchra 1977; Balzano 1983; Xu & De Zotti 
1989; Mazzarella et al. 1991) and interacting galaxies (Lar- 
son & Tinsley 1978; Kennicutt et al. 1987; Telesco et al. 
1988; Xu & Sulentic 1991) have on average significantly 
higher current star-formation rates than standard galax- 
ies like the Milky Way. Therefore we select these galaxies 
for our starburst galaxy sample. It includes those galaxies 
classified in the UGC catalog (Nilson 1973) as peculiars, 
double, trible, or multiple systems, and galaxies identi- 
fied by Markarian and collaborators as having ultraviolet 
excesses (see Mazzarella & Balzano 1986 for a review). 
Galaxies with active galactic nuclei (e.g. Seyfert galax- 
ies) are not considered. The region of low galactic latitude 
(|b| < 30°) is avoided in order to minimize the influence of 
Galactic extinction. Virgo galaxies listed in the VCC cata- 
log (Binggeli et al. 1985) are excluded, because in clusters 
galaxies may have additional radio emission mechanisms 
other than the star-formation alone (Gavazzi et al. 1991; 
V61k & Xu 1994). There are 133 objects in the starburst 
galaxy sample. 

The control sample includes spiral galaxies which are 
at high Galactic latitude (|b| > 40°), and from the mor- 
phological types Sa through Sd in the UGC catalog. 
The sample excludes Virgo cluster, Markarian, as well as 
Seyfert galaxies. Galaxies catagorized simply as 'S' are not 
included, because the probability that they are misclassi- 
fied early types is relatively high. There are 397 spirals in 
the control sample. 

For the majority of objects in both samples, fluxes at 
60/im and 100/im are taken from the Second Version IRAS 
Point Source Catalog. The exceptions are those extended 
(size > 5') and faint (feo^m < 0.6 Jy) sources listed in 
the catalog of IRAS Faint Source Survey (FSS) fluxes of 
CfA galaxies (Thuan & Sauvage, 1992), for which the data 
are taken from Thuan & Sauvage (1992). For 133 objects 
in the starburst sample, 113 are detected at both 60yum 
and lOO^m, corresponding to a detection rate of 85%. For 
the control sample, 329 are detected out of 397 objects, 
i.e. a detection rate of 83%. The integrated flux in the 
wavelength range 42.5 — 122.5 /im, Sfir , or its upper 
limit, is calculated from the IRAS fluxes (limits) accord- 
ing to Helou et al. (1988). Radio fluxes at 2.4 GHz are 
taken from the Arecibo survey data (Dressel & Condon 
1978). Fourty-six objects in the starburst sample and 92 
in the control sample are detected in radio, correspond- 



ing to detection rates of 35% and 23%, respectively. For 
an overwhelming majority of sources in our samples, the 
total "face-on" blue magnitude B^, corrected for Galactic 
extinction, inclination effect and for redshift, can be found 
in the Third Reference Catalog (RC3) (de Vaucouleurs et 
al. 1991). For a few exceptionals whose B^, are not avail- 
able, the blue magnitude is simply taken from the pho- 
tographic magnitude given in UGC and corrected for the 
Galactic extinction calculated using the formula given in 
RC3. Taking H = 75kms _1 Mpc _1 , the distances are cal- 
culated from the recession velocities after correcting for 
the Virgocentric flow (Aaronson et al. 1982). 



3. Statistical results 



Fig. la. The FIR and the radio luminostities for the starburst 
sample. Upper limits are indicated by arrows. 



In the statistical analysis carried out in this section, 
we use the so called 'survival technique' (Schmitt 1985; 
Feigelson & Nelson 1985) to exploit the information car- 
ried by the upper limits. 

We plot in Figs, la and lb the logarithms of the ra- 
dio intensities versus the logarithms of the FIR luminosi- 
ties for the starburst sample and for the control sample, 
respectively. The crosses are detected data, and the ar- 
rows denote the upper limits. Both samples show tight 
correlations between the FIR and the radio luminosities: 
the linear correlation coefficients are 0.90 and 0.88 for the 
starburst sample and the control sample, respectively. The 
slope of the log-log plot of the starburst sample is very 
close to unity: 0.98±0.01. On the other hand, the radio 
intensities increase slightly faster than the FIR luminosi- 
ties for the normal spiral galaxies in the control sample: 
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Fig. lb. The FIR and the radio luminostities for the control 
sample. Upper limits are indicated by arrows. 

the slope is equal to 1.09±0.01 in Fig. lb. In both plots, 
no galaxy is detected to seriously violate the correlation. 

Similar to the q parameter defined by Helou et al. 
(1985), we introduce a parameter 92.4GHz: 



Fig. 2. The logarithmic ratio of the FIR-to-radio luminostity, 
(/2.4GHz, and the indicator of the starburst strength, SB, defined 
in Eq. 2, for the starburst sample 



q 2 .4GHz = log[(S Fm /3.75 10 12 )/S„(2.4GHz)] 



(1) 



For the starburst sample, the mean of the g2.4GHz is 2.40 ± 
0.03, while for the control sample it is 2.52 ±0.02. The sta- 
tistical dispersion of the logarithm of the FIR/radio ratio, 
calculated using the survival technique, is a — 0.22 for 
the starburst sample, and a = 0.21 for the control sam- 
ple. Therefore, although the (?2.4GHz value of the starburst 
sample is slightly lower than the control sample, the dif- 
ference is within the dispersion. 

An important question in this study is: does the 
FIR/radio ratio depend on the starburst strength? In prin- 
ciple, the starburst strength can be defined either via 
the FIR-to-blue luminosity ratio or via the radio-to-blue 
luminosity ratio. Both reflect, approximately, the ratio 
between the star-formation rates averaged over the last 
~ 10 8 yr, s$, and the last ~3x 10 9 yr, sg, respectively 
(Xu et al. 1994). However, there is a trivial dependence 
(positive) of the FIR-to-radio ratio on the FIR-to-blue 
ratio, because the FIR luminosity is involved in both ra- 
tios. Similarly, a trivial dependence (negative) exists be- 
tween the FIR-to-radio ratio and the radio-to-blue ratio. 
In order to minimize this artificial effect, we introduce a 
starburst strength indicator SB invoking both ratios: 



SB =0.5 x [flog - (log 
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Fig. 3. <?2.4GHz and the logarithm of the FIR colour, log feo^m/ fioo^M 
for the starburst sample 

where (log (jf^j ) = -0.47 ± 0.02 and (log ( F Yg Hz ) } = 
— 15.70 ±0.05 are the means of the logarithms of the FIR- 
to-blue and of the radio-to-blue luminosity ratios of the 
control sample. A mean SB of 0.46±0.05 is found for the 
starburst sample, indicating that on average the sg-to-sg 
ratio is a factor of 2.9±0.4 higher than that of the nor- 
mal spirals in the control sample, consistent with previous 
studies (Kennicutt et al. 1987; Xu & Sulentic 1991). 
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Fig. 2 shows Q2.4GHz vs SB for 113 objects in the star- 
burst sample. The other twenty sources undetected either 
in the FIR or in the radio are dropped because no con- 
straint can be imposed on their (/2.4GHz values. No corre- 
lation is found in this diagram. The two galaxies with the 
highest SB values show very normal (72.4GHz- The varia- 
tion of the FIR-to-blue and the radio-to-blue ratios, so 
the SB, may be due to the variation of the extinction, in 
addition to the depence on the starburst strength. In or- 
der to verify the results obtained in Fig. 2, we plot in Fig. 
3 the dependence of (/2.4GHz on another - extinction- free - 
starburst indicator, namely the log/6o M m//ioopm, for the 
same galaxies as in Fig. 2. As pointed out by Xu & De 
Zotti (1989) and by Bothun et al. (1989), log fm^m/ fwo^m 
is tightly associated to the high mass star-formation rate 
of galaxies, and the galaxies with the strongest starbursts 
usually show the highest log f m)im / fion^-m values. Again, 
we do not detect any dependence of <?2.4GHz on this star- 
burst indicator. Therefore, we conclude that the strength 
of the starburst does not introduce any significant influ- 
ence on the FIR-to-radio ratio. 



4. The starburst model 

4-1. Star formation during a starburst 

In the model we describe the starburst galaxy as the sum 
of two components: a starburst population with a star for- 
mation rate (SFR) 'PsbW changing on short time-scales 
and a background population whose SFR is constant 
in time. The total SFR is the sum of both components, 
= *SB(t) + *b g - The SFR of the starburst population 
is parametrized in the following way: 



vj> (f ) = 1 *o exp*/ T «— for t < 
Sm ' ^* exp-*/ TSB fort>0 



(3) 



where T grow and tsb are the time-scales of the increase 
and the decrease of the SFR. The ratio between the maxi- 
mum SFR of the starburst population and the background 
population a = ^o/^bg is a measure of the starburst 
strength. An upper limit for a can be estimated from the 
total amount of molecular gas available. The SFR of the 
background population can be estimated by assuming that 
the total amount of gas present is used up in about 10 10 
yr. Then we can estimate an upper limit of a by: 



a = * /*bg < 



10 10 • x 

tsb 



(4) 



where x is the fraction of molecular gas available for the 
starburst. In a nuclear starburst, about 20 — 30% of the 
total amount of gas is in the nucleus and therefore avail- 
able for the starburst (Bernlohr 1992). The efficiency of 
the conversion of gas to stars can reach in starburst galax- 
ies 50% (e.g. Rieke et al. 1980). With these numbers we 
derive x = 0.1 — 0.15. Thus, e.g., for tsb 
upper limit for a is about 100. 



10 7 yr, an 



The Initial Mass Function (IMF), $(m), which deter- 
mines the mass distribution of the stars produced, is de- 
scribed by a power-law 



$(m)dm = m 7 dm 



(5) 



with an upper mass limit m u = 100 M Q and a lower mass 
limit m\ — 1 Mq (stars less massive than that are neither 
of importance for the FIR nor the radio emission). We 
adopt 7 = 2.7, in agreement with Scalo (1986) and the 
more recent results of Basu and Rana (1992). In addition, 
we considered the effect of changes in the slope of the IMF 
and m\. 



4-2. FIR emission 

The FIR luminosity is calculated as 



^FIR 



(*) = 



/' 

J m 



dm$(m)L(m)P(m)((m) 



r(m) 



dttt(t)(6)l 



where r(m) is the stellar main-sequenze life-time, L{m) 
the stellar luminosity (taken from Cox et al. 1988 and 
Cox et al. 1986), P(m) the fraction of light absorbed 
by the dust and C( TO ) the fraction of the absorbed light 
that is reemitted in the IRAS 40-120 /im band. Here, we 
parametrize the last two parameters in the same way as 
explained in detail in Lisenfeld et al. (1995): For mas- 
sive, ionizing stars (m > 20 Mq) we assume P(m) = 
0.6, £( m ) — 0.5. The radiation transfer of the non-ionizing 
emission of the intermediate massive stars (5 M Q < m < 
20 M Q ) and the old stars (1M < m < 5M ) is ap- 
proximated by an infinite parallel-slab geometry in which 
dust is homogenously mixed with the stars. We assume 
Ay = 1.5 for the starburst population (Bernlohr 1993) 
and Ay = 0.2 for the background population (Lisenfeld 
et al. 1995; Xu & Buat 1995). We make the further sim- 
plification that the radiation of the intermediate massive 
stars is absorbed according to the dust opacity in the UV 
(2000 A) and the radiation of the old stars according to 
the dust opacity in the blue (4400 A). 

Finally we split the FIR luminosity into the emission 
coming from H11 regions, L™^, and the diffuse emission 
ipiR' by assuming that stars more massive than 20 M 
are responsible for the dust heating in H11 regions, and 
less massive stars for the diffuse dust heating. 

4-3. Radio emission 

The radio emission consists of thermal bremsstrahlung 
and of synchrotron emission. The thermal radio emission 
is proportional to the total number of Lyman continuum 
photons ALyc emitted by the stars: 



PthermM = 1-08-10 



N- 



\ phots " 



Lyc 



GHz 



[W/Hz](7)l 
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(e.g. Condon 1992) where we have assumed that the frac- 
tion of Lyman continuum photons directly absorbed by 
the dust is / = 0.2 and the electron temperature in the 
Hll region T e = 7000 K. The total number of Lyman con- 
tinuum photons is calculated as: 



drn(j)(rn)hiy C (m) 



t—r(m) 



ditf(t) 



(8) 



where n-L, yc (m) is the production rate of Lyman continuum 
photons of a star of mass to, taken from Glisten & Mezger 
(1983). 

In order to calculate the synchrotron emission in a 
starburst a time-dependent equation has to be solved. 
The details are described in the Appendix. Here we will 
only outline the most important features. Several pro- 
cesses are taken into account: We assume that the CR 
electrons are accelerated in shocks of SN remnants. Then 
they propagate into the galactic disk and halo losing 
their energy radiativcly through inverse Compton and syn- 
chrotron losses. Since the radiation field in a starburst is 
very strong, the inverse Compton losses of the electrons 
are very high. This has the effect of shortening their radia- 
tive life-time, ti oss , and therefore of decreasing the tempo- 
ral shift between the FIR and the synchrotron radiation. 
(In a radiation field which is, for example, a factor of 10 
higher than in the Galaxy, i.e. about 10 eV cm~ 3 , the 
life-time is ti oss w 10 6 yr.) On the other hand, a high 
radiation field lowers the synchrotron emission itself, be- 
cause the number density of electrons at a given energy 
is lowered (see Appendix). Therefore, during a starburst 
- if the magnetic field stays constant - the synchrotron 
emission is a sensitve function of the radiation field and 
therefore of the starburst strength. In particular it is very 
low during the phase of most active star-formation. 

4-4- Starburst strength parameter SB 

In order to compare our model with the data displayed 
in Fig. 3, we compute in the model the starburst pa- 
rameter SB. The blue luminosity is obtained by assum- 
ing that the stellar emission is a blackbody spectrum. 
We obtain for the background, steady state population of 
the galaxy within our model: (log ^ L ™- ^ ) = —0.34 and 

(log ^ P2 ^ Hz ) ) = —15.48. The model predictions arc in 
rather good agreement with the mean values of the control 
sample. The small difference between model and data of 
(log (-^r 1 ) ), might be due to a slight overestimate of the 
average galactic extinction or due to the crude model as- 
sumptions (e.g. blackbody emission of stars). In the model 
-P2.4GHZ for the background population is calculated from 
the model result for Lfir and the mean value of (/2.4GHz 
of the control sample. Therefore the model prediction for 

(log ^ P2 ^° Hz ) ) is simply a combination of (log(-^Pp-)) 
(model), and <?2.4GHz (control sample). 



The value of SB depends very sensitively on the extinc- 
tion. The extinction in the starburst region is most likely 
higher than for the background population. We assumed 
Ay = 1.5, which is the typical internal extinction of the 
Hll regions in nearby galaxies (Kennicutt 1983), and at 
the same time consistent with the average Ay found for 
the Markarian galaxies by Xu & De Zotti (1989). 

In reality it seems plausible that some galaxies have an 
even higher extinction, explaining the large range in SB 
found for the starburst sample. The extinction mainly has 
an influence on SB, through Lb-, and hardly on the FIR 
emission (see Fig. 9) because even for the value of Ay 
assumed for the background population (Ay = 0.2) the 
galaxies are optically thick for the nonionizing UV light 
which is the main heating source for the dust (Xu 1990). 

5. Variation of q2.4GHz during a starburst 

In order to show the general features of the temporal 
evolution of (/2.4GHz during a starburst, we display in 
Fig. 4 the variation of g2.4GHz during a starburst with 
a timescale tsb = 10 7 yr and r grow = for different star- 
burst strengths, and in the Figs. 5 and 6 separately the 
temporal evolution of the FIR and the radio emission. 



Fig. 4. The ratio between the FIR and the radio emission at 2.4 
GHz, (/2.4GHz (i), is shown as a function of time for tsb = 10 7 
yr, Tg row = 0, B =const, a = ^o/^bg = 100 (full line), a — 50 
(dashed line), and a = 25 (dotted line). The dash-dotted line 
shows <?2.4GHz for a steady state galaxy. 

The temporal variation of (72.4GHz can be split into 3 
phases: In the early phase (phase I; here: until w 50 Myr) 
(?2.4GHz(i) is higher than the steady state value at this fre- 
quency. In this phase the emissions from H 11 regions dom- 
inate the total FIR and radio emission (see Figs. 5 and 6). 
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Fig. 5. For a starburst with a — 100 and otherwise the same 
parameters as in Fig. 4 the FIR emission (full line) is shown as 
a function of time. The FIR emission from H II regions (dashed 
line) and the FIR emission from diffuse dust (dotted line) are 
shown separately. The dash-dotted line describes the FIR emis- 
sion from the background galaxy. 



With respect to the FIR emission, this is due to the high 
relative abundance of massive stars. For the synchrotron 
emission other reasons are more important: First of all, 
it takes serveral 10 6 yr until the first SN explosions oc- 
cur. More important, however, is the fact that during the 
first serveral 10 7 yr the synchrotron emission is strongly 
suppressed (see Fig. 6) due to the strong inverse Comp- 
ton losses caused by a high radiation field. It is worth to 
note that only because of the presence of the thermal ra- 
dio emission and because the two ratios (log(L™^/P t herm) 
and log(Lp' n ^/P syn ) are relatively similiar no major differ- 
ence of <;2.4GHz with respect to a steady state galaxy occurs 
in this phase. 

In the next phase (phase II; here: between about 50 — 
100 Myr), (/2.4GHz can be even lower than in a steady state. 
This is due to the time delay between the SN explosions 
and the FIR emission: it leads to a maximum of the radio 
emission at a time when the FIR emission is already in its 
declining phase. 

In the last phase (phase III; here: after about 100 Myr) 
the opposite temporal effect occurs: The synchrotron ra- 
diation has already decreased to the steady state value - 
or even below because the inverse Compton losses are still 
higher than before the starburst - but the FIR emission is 
still high due to dust heating by low-mass (m = 1 — 3 M Q ) 
stars. Therefore during this phase g2.4GHz is higher than 
the steady state value and depends sensitively on a = 
^o/^bg- In this phase the galaxy is already completely in 
the post-starburst phase. 



Fig. 6. For a = 100 and otherwise the same parameters as in 
Fig. 4 the radio emission during the starburst is shown (full 
line), devided into thermal bremsstrahlung (dashed line) and 
synchrotron emission (dotted line). The dash-dotted line de- 
scribes the radio emission from the background galaxy. 



Fig. 7. The model result for the same parameters as in Fig. 4 
(with a = 100) are compared to the data. The full line shows 
the model results for SB and </2.4GHz for the time interval be- 
tween t = and t = 10 9 yr with the arrow indicating the 
direction of increasing time. 



The deviations of g2.4GHz from the steady state value in 
the various phases are therefore either due to the different 
time-scales of the FIR and the radio emission or due to the 
high inverse Compton losses experienced by the electrons. 
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In Fig. 7 the model results are compared to the data. 
It can be seen that in phase I the model predictions for 
92.4GHz lie substantially above the data. In the following 
we will try to find parameters that reduce the temporal 
variation of Q2.4GHz to a minimum and fit the data. 



6. Constraints on starburst parameters by 
comparison with the observations 

6.1. Young starburst (phase I): Fast increase of the 
magnetic field? 

The main reason for the high model value of g2.4GHz dur- 
ing the early phase is a lack of synchrotron radiation due 
to high inverse Compton losses. The radio emission during 
the first ~ 2 x 10 7 yr is therefore - if the magnetic field 
remains constant - practically completely thermal (Fig. 
6). This has not been observed, however: Van den Driel et 
al. (1991) determined the thermal radio fraction of a sam- 
ple of starburst galaxies and derived only 40% at 5 GHz 
(corresponding to about 30% at 2.4 GHz). A possibly very 
high thermal fraction should furthermore be noticeable in 
a flat radio spectrum (P rac j(^) oc v~ QA ) which has not 
been observed either. 

The suppression of the synchrotron emission will be 
substantially diminished if the magnetic field increases 
rapidly already in phase I. The behaviour of the mag- 
netic field in a starburst has not been well studied yet, but 
it seems plausible that the field increases. Ko & Parker 
(1989) proposed a model in which the increase of inter- 
stellar turbulence due to a high SFR activates a turbulent 
galactic dynamo. An alternative process to increase the 
magnetic field in a starburst could be the compression by 
SN shocks of the ionized gas in which the field is frozen. A 
third field enhancement process could be the systematic 
dynamo effect due to a commencing galactic wind as we 
may be observing in M 82 or even at the centre of our 
Galaxy. 

In each of the three scenarios the increase of the mag- 
netic field is ultimately caused by the star formation. Af- 
ter the end of the starburst, the magnetic field enhance- 
ment stops and the field dissipates slowly. Therefore we 
will tentatively parametrize the temporal variation of the 
magnetic field as: 



B{tf 



oc £fir(*) 
= constant 



if 



dL-B 



dt 

otherwise 



> 



(9) 



In Figs. 8 and 9 the results are shown. It can be seen 
that the magnetic field increase lowers substantially the 
value of 92.4GHz in phase I and allows to fit the data. 

We would like to stress that the results presented do 
not crucially depend on the exact parametrization of the 
magnetic field increase during phase I. The important 
point is a fast increase of the magnetic field. The time- 
scale of this increase has to be rather short (« 10 7 ) be- 
cause a large magnetic field strength is required already 



Fig. 8. Q2.4GHz is shown for r SB = r grow = 10 7 yr, a = 100 



(dashed line), tsb = Tgrow = 3 x 10 7 yr, a = 50 (full line) 

3 x 10 7 yr, 



(both with mi = IMq) and for tsb = T gr0 w 
a — 50, mi = 3 M (dotted line). An increase of the magnetic 
field according to Eq. (9) has been taken into account in all 
three cases. 



Fig. 9. In the same way as in Fig. 7 the model results and 
the data are compared. The parameters for the model are: 
Tgrow = tsb = 3 x 10 7 yr, a = 50, mi = 1M , Ay = 1.5 (full 
line). In the dashed line the result for Av = 5 is shown in order 
to illustrate the influence of the extinction on the parameter 
SB. 



in an early stage of the starburst. The time-scale for the 
dynamo activity has been estimated to be between 10 8 
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(Ko & Parker 1989) and > 10 9 yr (Field 1993). Therefore 
a dynamo may be a too slow process for a starburst. How- 
ever the alternative processes, such as outward "combing" 
of the field by a mass outflow or of field compression by 
SN shocks, would occur on the time-scale of stellar wind 
bubble evolution and SN activity. At least the wind ac- 
tivity is directly coupled to the FIR enhancement with 
delays of the order of a million years only. An example for 
this might be the nearby starburst galaxy M 82. It lies on 
the FIR/radio correlation (e.g. Wunderlich et al. 1987). 
Bcrnlohr (1992) estimated that the starburst in M 82 has 
been active for 1 — 5 x 10 7 yr, i.e. M 82 is in our phase I. 
The fast galactic outflow from the starburst nucleus and 
the corresponding convective electron transport is consis- 
tent with the observed synchrotron emission and the SN 
rate (Volk et al. 1989). 

Parker (1992) has argued that galactic dynamos can 
be driven by the pressure of cosmic rays and the sub- 
sequent reconnection of azimuthal field lines to produce 
poloidal fields. This process is competitive to our assumed 
extension of magnetic field lines to "infinity" by a wind. 
The point here is that also the above version of a galac- 
tic dynamo can more or less adiabatically follow the time 
variations in the production of cosmic rays, presumably in 
SNRs, as our data indeed suggest. 



6.2. Evolved starburst (phase II): Time-scale 

If the time-scale of the variation of the star formation rate 
is very short, i.e. roughly shorter than the time-scale for 
SN-explosions (~ 10 7 yr, i.e. the life-time of SN progen- 
itors), the maximum of the attendant synchrotron emis- 
sion occurs at a time when the FIR emission, which has 
about the same time-scale as the starburst, has already 
decreased significantly. This results in a very low model 
value of <72.4GHz (less than 3 a below the steady state 
value). Therefore, in order to fit the data satisfactorily 
during this phase, T grow + tsb has to be at least compara- 
ble to the SN time-scale, i.e. about 10 7 yr. Such time-scales 
are short enough to allow for repetitive bursts of star for- 
mation on time-scales of 10 8 yr as discussed by Kriigel 
& Tutukow (1993). Our model results provided good fits 
to the data for time-scales longer than about 3 x 10 7 yr. 
These numbers are in agreement with other authors who 
found starburst time-scales between 2 x 10 7 and some 10 8 
yr (Rieke et al. 1980; van den Broek et al. 1991; Bcrnlohr 
1993). We furthermore inferred a minimum time-scale for 
the increase of the SFR of T grow ^3x 10 7 yr. This slow 
increase was necessary in order to provide a large enough 
magnetic field strength and SN rate in the early starburst 
phase. 



6.3. Post- starburst (phase III): A higher low-mass cut-off 
of the IMF? 

We also tested the effect of changes in the IMF, both of 
the slope k and the low-mass cut-off m\. These two pa- 
rameters are suggested by many authors to be different in 
a starburst in the sense that the star formation is likely 
to be biased towards massive stars (e.g. Rieke et al. 1980; 
Rieke et al. 1993; Bernlohr 1992; Wright et al. 1988). For 
our model, the most noticeable influence of the IMF oc- 
curs in the post-starburst phase (III). When the canonical 
IMF (Eq. (5)) is assumed for the starburst, the enhance- 
ment in the FIR emission lasts significantly longer than 
the radio emission because of the dust-heating due to the 
A stars (~ 2M©) produced in the burst. This is shown 
by a relatively high model value of (72.4GHz in this phase 
(Figs. 4 and 8). A better agreement with the data can be 
obtained if a higher low-mass cut-off, i.e. m\ ~ 3 M Q (see 
Fig. 8), is adopted. A flatter slope will have a similar ef- 
fect. Hence, our result seems also in favour of a top-heavy 
IMF for starbursts. 



7. Summary and Conclusions 

We have analysed and compared the FIR and radio data 
of a sample of interacting/Markarian galaxies (starburst 
sample) and a sample of normal galaxies (control sam- 
ple) in order to test for the influence of the star formation 
activity on the FIR-to-radio ratio. A model for the FIR 
and the radio emission in a starburst was used to inter- 
prctc the observational results and to derive constraints 
on starburst parameters. The following main conclusions 
could be drawn: 

The mean FIR-to-radio ratio, (72.4GHz, showed no sig- 
nificant difference between the starburst sample and 
the control sample and is independent of the starburst 
strength. 

The magnetic field is required to increase strongly on 
short time-scales (« 10 7 yrs) in order to maintain a 
sufficiently constant (72.4GHz m the first « 2 x 10 7 yr 
of a starburst. Otherwise the inverse Compton losses 
during this phase would lead to a very low synchrotron 
emission, resulting in a value of (72.4GHz significantly 
higher than the observed one. 

The time-scale of the variation of the star-formation 
rate has to be longer than some 10 7 yr. Otherwise the 
different time-scales of the FIR and the radio emission 
produce large fluctuations in (72.4GHz- 
Our results favour a top-heavy IMF for starbursts (i.e. 
either a higher lower mass cut-off of m\ ~ 3M or a 
flatter slope of the IMF) because this reduces the con- 
tribution of A stars to the dust heating, thus short- 
ens the time-scale of the FIR emission and therefore 
improves the agreement with the data for the post- 
starburst (phase III). 
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In order to calculate f(E,t) from Eq. (A2), one has 
to know ^sn(*) and b(t), the latter being proportional 
to t/ rac i(£) + Usit). The supernova rate ^sn(*) is calcu- 
lated directly from the SFR given in Eqs. (3) and (5). The 
temporal evolution of the energy density of the radiation 
field, U Ta d(t), is proportional to the bolometric luminosity 
^boi(i) of the galaxy if the volume emitting the radition 
does not change significantly during the starburst: 



Appendix 
burst 



Synchrotron emission during a star- 



U Ta .d(t) = L bo i(t) 



r/ bg 

U rad 
7-bg 
^bol 



(AA) 



Assuming that the time-scale of the radiative energy losses 
is much shorter than the time-scale of diffusion, which 
is certainly a valid assumption for a starburst where the 
inverse Compton losses of the electrons are very large, we 
can neglect diffusion and describe the temporal evolution 
of the electron particle density f(t,E) by the following 
equation: 



df(t,E) 

at 



E 



qSNV^(t) + -^{hU ■ I; 2 J U.E) }(A1 )| 



where (E/mc 2 )^ x qsN^SN(t) is the source spectrum of the 
relativistic electrons with i/sn(*) as the supernova rate and 
qsN as the number of relativistic electrons produced per 
supernova and per energy interval. The power law index of 
the source spectrum is taken as x — 2.2 (Volk et al. 1988). 
The energy losses b{t) oc C\(JJ Ta d(t) + Us(t)) describe the 
inverse Compton and synchrotron losses which are propor- 
tional to the energy density of the radiation field C/ ra d and 
the energy density of the magnetic field Ub, respectively. 

The solution of Eq. (Al) for the starburst population 
is given by: 



We take the radiation energy density of the steady state 
galaxy before the starburst U ' J*, — leV/cm 3 , according to 
the value in our Galaxy. i£oi> its bolometric luminosity, is 
derived from the model, assuming a constant SFR. 

The temporal evolution of the magnetic field is parametrized| 
according to Eq. (9). The energy density of the magnetic 
field in a steady state is taken as the Galactic value: 
ifg S = leV/cm 3 which corresponds to a magnetic field 
)lf 5 bg = 6.3^G. 

' Knowing the electron particle density f(E,t), we can 
calculate the synchrotron radiation by: 



dE 

P syn (v,t)dv = f(E(v),t) — 



syn 



dE 

cU7 dl/ 



(A5) 



with (dE/dt)\ syn as the energy loss by synchrotron radia- 
tion and dE/dv given through: 



E \v G {B) 

2 



mc 



0.45 



(A6) 



/SB (£,*) = 



E 



■ <ZSN 



L 



max(t s tart,t — Tl ss(E,t)) 



dt' V ^{t') [l-E 



j b{t")dt" 



x-2 



(A 



Here, ri oss (.E, t) is the time-scale of the radiation losses of 
the electrons during the starburst which is defined by the 
relation: 



(Longair 1994) with vq being the gyrofrcquency. This cor- 
iponds to the simplification that an electron of energy 
lemits the whole synchrotron radiation at the maximum 
^quency of the synchrotron spectrum. 

In order to find the absolute normalization of the syn- 
Irotron radiation we scale the synchrotron radiation such 
that the FIR/radio ratio (?2.4GHz for the background popu- 
lation is equal to the observed average value of the control 
sample. 



E 



t-Tioss(E,t) 



b(t')dt' = 1. 



(A3) 

References 



The boundary condition is that at the beginning of the 
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